Effect of perturbations on information transfer in spin chains 
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Spin chains have been proposed as a reliable and convenient way of transferring information and entanglement 
in a quantum computational context. Nonetheless, it has to be expected that any physical implementation of 
these systems will be subject to several perturbative factors which could potentially diminish the transfer quality. 
In this paper, we investigate a number of possible fabrication defects in the spin chains themselves as well as the 
effect of non-synchronous or imperfect input operations, with a focus on the case of multiple excitation/qubit 
transfer. We consider both entangled and unentangled states, and in particular the transfer of an entangled pair 
of adjacent spins at one end of a chain under the miiToring rule and also the creation of entanglement resulting 
from injection at both end spins. 

PACS numbers: 03.67.-a, 03.67.Lx, 75.10.Pq 
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I. INTRODUCTION 

A crucial ingredient in quantum information processing 
based on solid state systems is the transfer of quantum in- 
formation. Assuming that there are quantum registers for 
computing and storing information, the ability to transfer this 
information reliably and efficiently from one register to an- 
other is vital for the construction of larger, distributed and 
networked systems. A solution to this challenge has been pro- 
posed through the use of spin chains ||11|31- The mathematical 
framework underpinning spin chains can be applied to vari- 
ous physical devices; these could be made of any components 
whose states can be mapped onto spin i particles interacting 
with their neighbors. Electrons or excitons trapped in nanos- 
tructures form explicit examples 13|-|5|], as do nanometer scale 
magnetic particles la] or a string of fuUerenes 131. Another 
representation is the encoding into a soliton-like packet of ex- 
citations ill]. 

Within spin chains, a single-site excitation is defined as an 
"up" spin in a system that is otherwise prepared to have all 
spins "down". A discussion about unmodulated spin chains 
has been given in 13, uM whereas in iflHl the couplings were 
chosen to be unequal. There has also been research on wire- 
like chains with controlled coupling strength at either end lll2ll 
and transfer through parallel spin chains lllSll . to name but a 
few closely related areas. Here we only consider linear spin 
chains whose coupling strength Ji^i+i between two neighbor- 
ing sites i and i + 1 has been pre-engineered to ensure perfect 
state transfer (PST) along the chain ll5L [l4ll . For a chain of 
length N with characteristic coupling constant Jq, the PST 
coupling strength sequence is defined as Ill4ll 
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For devices based on excitons in self-assembled quantum 
dots, Jo is mainly governed by Forster coupling |3], which 
in turn depends on the distance between the dots as well as 
the overlap between the electron and hole wavefunctions in 
each dot. In gate-defined quantum dots, however, Jg will de- 
pend on tunnelling and thus on parameters such as the width 
and height of the barriers which separate the different dots, as 
well as on the overlap of electronic wavefunctions centered in 
different dots. For chains of fullerenes or actual atoms Jq will 
represent some "hopping" parameter describing the propen- 
sity of the excitation to transfer from one site to the other. The 
natural dynamics of a spin chain can then be described by a 
time independent Hamiltonian as follows 



N 



N-1 



n = ^e,|l)(l|, + J2 -^M+iIllXOl. ® |0)(l|m + 

|0)(l|,®|l)(0|,+i]. (2) 

In a perfect system (to which perturbations will then be ap- 
plied) we will assume the single excitation energies e^ to be 
independent of the site i, and therefore only concentrate on 
the second term of Eq. (|2]l. In some physical systems such 
as quantum dot strings, e^ could naturally differ according to 
position, but may be tuned to be the same at all sites via ap- 
plication of local fields yl]. The fidelity F, corresponding to 
mapping the initial state |V'm) over a time t into the desired 
state iV'/in) by means of the chain natural dynamics, is given 
by 



F = \{iJf^n\e-'^'^''\^^n^)\' 



(3) 



and PST is realized when the evolution is arranged to achieve 
F = I. We use the fidelity of state vectors to determine 
the transfer quality of information for unentangled states, as 
detailed for example in yj]. For entangled states, we mea- 
sure instead the entanglement of formation (EoF) as defined 
inRef. UJ. 

The time evolution of a system is dependent on its characteris- 
tic coupling constant Jq. In particular, the time scale for PST 



from one end of a chain to the other, also known as mirroring 
time, is tjv/ = 7rfi,/2 Jq so that the periodicity of the system 
evolution is given by tg = ttH/Jq. As the Hamiltonian ^ 
preserves the excitation number, the evolution of the initial 
state will remain within the original excitation subspace. 



II. EFFECT OF FABRICATION DEFECTS 

We will now consider the influence of general fabrication 
defects on linear spin chains with multiple excitations. 

(a) Random noise 

We model the effect of fabrication errors (random, but fixed in 
time) for the energies and couplings in the system by adding 
to all non-zero entries in the Hamiltonian matrix a random en- 
ergy rjdi^mJo for 1 < l,m < number of basis states. The 
scale is fixed by 77 which we set to 0. 1 and for each I < m 
the different random number d/ ,„ is generated with a flat dis- 
tribution between zero and unity. For the other side of the 
diagonal with m < I, di^„i = dmj, preserving the hermitic- 
ity of the Hamiltonian. This method of including fabrication 
defects means that we could observe effects of a reasonable 
magnitude although clearly other distributions could also be 
modeled; for specific tests, the weight of the noise would have 
to be determined according to the individual experiment being 
simulated. 

(b) Site-dependent "single-particle" energies 

As a further possible fabrication defect, we consider the ef- 
fect of the first term of (|2]i that we previously dismissed under 
ideal conditions 
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Hi may represent external perturbations, such as local mag- 
netic fields, or additional single site fabrication imperfections. 
We thus assume here that e^ is not independent of the site / 
any more. 

(c) Excitation-excitation interactions 
In spin chains with multiple excitations, we also consider a 
perturbation term 
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which represents the interaction between excitations in 
nearby sites. For example, this may correspond to a biexci- 
tonic interaction in quantum dot-based chains 11 ml 1711 . 
(d) Next-nearest neighbor interactions 
Finally, we also investigate the effect of unwanted longer 
range interactions, which could be an issue when consider- 
ing pseudospins based on charge degrees of freedom. For this 
we add to ^ the perturbative term 

i?3= E ^M+2[|l)(0|««'|0)(l|.+2 + 

|0)(l|,®|l)(0|,+2]. (6) 



The expression for Ji.i+2 will depend on the type of inter- 
action between spin chain sites. Here we explicitly consider 
three cases. The first and more general approximates the next- 
nearest neighbor interaction as proportional to the average of 
the related interactions between nearest neighbor sites. 



Ji,i+2 — ^{Ji,i+1 + «/? + !, 1+2)72, 



(7) 



with the parameter A defining the strength of the interaction. 
This expression simulates the original coupling modulation of 
the chain. 

Secondly we explicitly consider dipole-dipole interactions, 
which are relevant, e.g. to chains of quantum dots with exciton 
qubits and Forster coupling lllSll . In this case the coupling 
between sites scales as 1/R^, with R the distance between the 
two sites considered. For roughly equidistant sites, we then 
expect the next-nearest neighbor couplings to be about a tenth 
of the nearest neighbor couplings. By using this and Eq. ([B 
we obtain 



j: 



dip 
iA-\-2 



J( 



{[^(A^ - t)]-^ + [{t + 1){N - I - I)]--- }' 



(8) 

Finally we consider the case of coupling due to tunnel ling , 
relevant e.g. to graphene quantum dots with spin qubits 111911 . 
Here the coupling scales as At^/U, with U the on-site 
Coulomb energy and t oc iie^^l'^l the tunnelling (hopping) 
parameter, with k the 'forbidden' momentum in the bar- 
rier III9I1 . In this case there is no explicit expression for 
the next-nearest neighbor couplings in terms of {Ji.i+i}, but 
{JII.+2} '-^'1 ^^ determined numerically by using the expres- 
sion for t and Eq. ([T]i- We note that we expect { Jj*™2} '^o ^^ 
very small, as the interaction decays exponentially with the 
distance. 

If we consider chains of quantum dots with exciton 
qubits, we can assume J max ~ ImeV, where J„iax = 
maXi{Ji,i+i}, and Rmin ~ lOnm. Results in Fig. [T]then 
show that this is in very good agreement with considering Eq. 
(|7]i with A = 0.12, a value that one would expect from dipole- 
dipole interaction and basically equidistant sites (see discus- 
sion above). In this case the effect of next-nearest neighbor 
interaction is extremely detrimental to the system as even the 
fidelity of the first state transfer at t = Q.5ts is reduced by 
almost 50%. 

By contrast, the tunnelling mechanism case leads to very 
different results. Using the parameters in 111911 for graphene 
dots and spin qubits, we obtain the results shown in the inset 
of Fig. [T] after 3ts the system has lost less than 2% of its 
fidelity. These results are in very good agreement with Eq. (JTl 
with A = 0.0001 as can be seen in the inset of Fig.l, where 
the plots are virtually indistinguishable. It is very encouraging 
that in this case realistic parameters point to such small values 
of A and thus generate very high fideUty transfer 

In the following we will use the expression in Eq. O to 
further discuss the effects of H3. As values of A smaller than 
0.01 have a minor detrimental effect on the system, we will 
from here onwards focus on the range 0.01 < A < 0.1, where 
the upper bound may be of interest to some experimental im- 
plementations, e.g. as discussed in the dipole-dipole interac- 
tion case. 
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FIG. 1 : Influence of H3 on the example of an 8-spin chain with initial 
input state jV'in) = j 11000000), fidelity vs. rescaled time t/ts, 
according to both Eq. {Tjl and Eq. ([8} (dipole-dipole coupling). The 
peak at t = 0.5fg is approximately 0.52. As the model of Eq. Q 
matches the data derived from Eq. {Sj extremely well, the lines of 
the respective plots are nearly indistinguishable. 
Inset: As for main panel but for graphene quantum dots and spin 
qubits, with tunnelling coupling (thin lines) and coupling according 
to Eq. (|7](, A = 0.0001 (bold lines). Again results from the two 
models for the coupling constants are almost indistinguishable from 
each other. 



We will now consider the effect of the fabrication defects 
(a) to (d) first on the transfer of factorisable states, i.e. unen- 
tangled chains, and then on entanglement creation and entan- 
glement transfer along a spin chain. For unentangled states, 
we will consider a 6-spin chain for all investigations of fab- 
rication defects while in the case of entanglement, we also 
consider an 8-spin chain. As we will explicitly see later for 
Hi, the influence of fabrication defects does depend on the 
chain length, but not on the parity of the chain. 



A. Transfer of unentangled states 



FIG. 2: Influence of random noise (see (a) in text) for r; = 0.1 on 
a 6-spin chain with two excitations, fidelity vs. rescaled time t/ts- 
States other than the input state and its twin state are not shown. The 
first peak at f = 0.5ts = tm is 0.9975. 
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FIG. 3: Effect of perturbation H2 (Eq. l|5j) on a 6-spin chain, fidelity 
vs. rescaled time t/ts- 

(a) 7 = 0.05: the first peak at t = 0.5ts = tm is 0.9988. 

(b) 7 = 0.1: the first peak at t = 0.5is = tM is 0.9954. 



The device we consider in this paper is a linear spin chain 
with couplings fixed such that the conditions for PST are sat- 
isfied. One of the properties of these devices, which we make 
heavy use of, is the mirroring rule 1I20I ulll . The mirroring 
rule is such that any state injected into a linear spin chain, 
subject to the constraints of Eq. ([T]i and site-independent ci, 
evolves into its "mirror state", where the symmetry center of 
a chain is the middle point in chains with N even, or middle 
site {N + l)/2 in chains with N odd. This mirroring property 
is independent of the number of excitations a state comprises 
and also of the length of the chain. Furthermore the mirror- 
ing rule holds for states spread across the whole chain: all 
excited sites are mirrored across to their "twin" with respect 
to the chain center of symmetry. We call this mirrored state 
the "twin state". When investigating the quality of a device 
with fabrication defects of any sort, we are primarily inter- 
ested in the effect on the "twin state" att = tM but it is also 
important to consider the effects on the next few periods as 



some quantum information protocols may require periodicity 
of their systems PZ]. 

As can be seen in Fig. |2] the effect of random noise leads 
to a continuous and definite decline in the transfer fidelity of 
both the input state |110000) and its "twin state" |000011) at 
the mirroring times. Naturally, this means an increased prob- 
ability of the occurrence of other possible states (not shown 
on the graph) but this happens on a relatively unpredictable 
basis, with no one state ever becoming and remaining partic- 
ularly prominent. 

In comparison, when looking at the evolution over a few 
periods, even a relatively large value of 7 = 0.05 in Eq. ^ 
has less effect on the 6-spin chain, as shown in Fig. [3] (a). 
Similarly, the on-site energies represented by Eq. (JUi also lead 
to an unrecoverable decay in the state transmission fidelity, 
with the excitations being ultimately entirely spread out along 
the chain (not shown). A combination of any of these three 
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FIG. 4: (Color online) Revival fidelity of the state ]?/;,„) = 1 111000) 
in a linear 6-spin chain, measured on the first revival peak at time 
t = ts, vs 7 and e . 
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FIG. 5: Influence of H3 (Eq.©) with A = 0.01: fidelity of the state 
l^in) ~ |110000) in a 6-spin chain vs. rescaled time t/ts- The first 
peak at t = 0.5is = ti\/ is 0.9966. 



perturbation factors simply accelerates the decay trend of the 
desired states. 

However, if we simply analyze the first revival peak at t ^ ts, 
we see that for varying 7 and e, even when perturbing the 
system by as much as 20% of Jo, the system suffers less than 
a 10% loss in fidelity (Fig. |4]i. Here, to represent the fact that 
Ei is site-dependent, we used a value e^ = eJodi for all sites 
/, with < di < 1 a random number generated from a flat 
distribution for each i. An average over 200 realizations for 
every value of e was then taken. 

Finally, we also note that next-nearest neighbor interactions 
perturb the system similarly to noise. This can be seen in Fig. 
|5]where even for a reasonable value of A = 0.01, the transfer 
peaks of both the initial state of a 6-spin chain and its mirror 
twin quickly decay. In order to achieve the same fidelity loss 
at i = 4.5is as for A = 0.01 (Fig. |5]), we have to consider a 
value of 7 one order of magnitude bigger, as is shown in Fig. 
|3](b). 

In our simulations we have kept the value of Jmax constant 




FIG. 6: (a) Fidelity of initial vector |110 ■ • ■ 0) at i = is vs. chain 
length A'^ for three values of A, as labelled. Fits are according to 
e-^^'Mi with Ao = 0.21. 

(b) Average fidelity of initial vector jllO ■ ■ ■ 0) at i = fs vs. chain 
length A^ for three values of e, as labelled. Fits are according to 
e-'^'''/'l witheo = 1.12. 



as N is varied: this models the physical constraint that in any 
realistic system the coupling strength is capped by a maxi- 
mum characteristic value. As a result Jq ~ 2Jmax/N ( Jq = 
'2'Jmax/N ^\ — l/iV^) for even (odd) chains. To avoid this 
implicit dependence on N, we have here set ei = eJmaxdi 
(where di is a random number from a uniform distribution 
within and 1) and averaged every point on the graph from 
200 random realizations. As we see in part (a) of Fig. |6] the 
effect of Hs becomes very detrimental to the system even for 
the relatively small value of A = 0.05, although the loss of 
fidelity for A = 0.01 is very small even for long chains. On 
the other hand, the effect of on-site energies may be tolerable 
for values of e up to 0.1, where long chains of A^ = 15 suffer 
less than a 20% loss in fidelity. 

A simple estimate of the effect of errors in the energy level 
spectrum u4M suggests an overall error, or loss in fidelity for 
PST that scales as an exponential decay in N with Gaussian 
dependence on the characteristic noise parameters. We com- 
pare this with numerical results in Fig. |6] where the loss 
of fidelity due to Hi scales as e~^^^'^'> and the loss due to 
H3 scales as e"^^^'^^ with increasing chain length N, where 
/(A) = A^/Aq and /(e) = e^/eg so that Aq and eo charac- 
terize the impact of the noise. The comparison in Fig. |6]shows 
that this simple analytical form can reproduce the numerical 
results to a high degree of accuracy: the loss of fidelity scales 
indeed as an exponential decay with Gaussian damping in the 
noise parameters. 

We conclude from these results that the transfer of unentan- 
gled states across the device is very robust against the pertur- 
bations Hi and H2, and less so for H^. The effect of noise 
in the system, as we have implemented it, is slightly more 
noticeable at the transfer time t]\f but still allows for excel- 
lent transfer at a loss rate of just over 10% over the course of 



4 periods. With regards to longer term periodicity, it is the 
next-nearest neighbor interaction term H^ which perturbs the 
system most for the values of e^, 7 and A shown. 



B. Transfer and creation of entangled states 

One of the most outstanding properties of spin chains is 
their ability to not just transfer reliably factorisable states, but 
also to transfer information encoded as entangled states. This 
is again based on the mirroring rule 1I20II21I1 and was first men- 
tioned and discussed under various aspects in Refs. lll4l l23l - 
I28II . Entanglement is one of the key resources in quantum 
computing, and is crucial to some quantum cryptography pro- 
tocols, and to quantum teleportation. Being able to reliably 
transfer entanglement from one place to another is therefore 
a core interest that the device we are analyzing should be 
able to respond to. This property follows from the fact that 
the set of states |(/){i}_„) = [ii,i2, • • • , Jiv) - *j == 0, 1 to 
which the mirroring rule applies, is a basis set so that any 
state, and in particular entangled states, can be written as the 
superposition |*) = E{,;},„ C{j},nl'/'{j},n)- Here n is the 
total number of excitations in the state, and {1} the ensem- 
ble of indices ij different from zero. The above relation im- 
plies that after t^i has passed the following "twin state" is 
reached: |*toi„) = !]{»},, i C{,},„exp{-ii?iM}|<?!>{t},«) = 
^{i}.n C{i},ri|0{i},ri;t«,in), SO that in particular any entangled 
state is transferred into its mirror entangled state. 

When considering spin chains with (i) an entangled initial 
state or (ii) an initial state that leads to entanglement, it is more 
useful to observe the evolution of the EoF in the system. We 
use chains with an initial Bell state on spins 1 and 2 to repre- 
sent scenario (i), as for example j^/im} = 4=(|1000) + |0100)). 
Accordingly, we monitor the EoF in the reduced density oper- 
ator of spins (1,2), tracing out the rest of the chain. To mon- 
itor the "twin state" entanglement we calculate the EoF for 
the reduced density operator of spins {N — 1,N). Scenario 
(ii) is different in that the initial state of the chain is not en- 
tangled, but will lead to entanglement through natural dynam- 
ics. As an example of this, we use a linear chain with input 
1+) = ^(|0) + |1)) on both spins 1 and A^. This is equivalent 

to an initial state |V'm) = lOiO^v + liOw + OiIat + IiIat) ® 
O2O3 • • • Oat-i), where the subscripts designate the spin site. 
A system set up in this way will then show a maximally en- 
tangled state in spins 1 and N at time t = Im 11221 I29l43lll . 

The effect of noise on entangled chains shows a similar 
trend to that of unentangled chains. However, we notice that 
the loss of EoF in Fig. |2]which represents case (ii) is nearly 
10% bigger than the loss of EoF in Fig. [8] (case (i)) over the 
course of the shown 7 periods. As a comparison, the unentan- 
gled state of Fig. [2]loses a similar amount of fidelity over the 
same amount of time as the chain in Fig. |7] 

There is no effect from perturbation H2, (Eq. (|5]i), in case 
(i), as there is only one excitation in the system in both am- 
plitudes. Obviously this would not be so for a Bell state with 
a doubly excited amplitude. Similarly, the effect of H2 in 
case (ii) is restricted to the two excitation subspace populated 
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FIG. 7: Influence of random noise (see (a) in text) for r) = 0.1 on 
the entanglement of formation between the end spins of a linear 8- 
spin chain with input state I VJin) = l/2(|0i08) + IliOg) + lOilg) + 
Ills)) ® IO2 ■ • ■ O7) vs. rescaled time t/ts- The first peak at t = 
0.5ts = tM is 0.9983. 
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FIG. 8: Influence of random noise (see (a) in text) for r; = 0.1 on 
the entanglement of formation between the end spins of an 8-spin 
chain with input state \ipi,^) = l/^/2(| 10000000) + JOIOOOOOO)) 
vs. rescaled time t/ts- The first peak at i = 0.5is = Im is 0.9993. 



by the evolution of the jlilw) <8i IO2 ■ • ■ Ojv-i) component of 
\'4)in) only and is thus not very prominent, as is shown in Fig. 

121 

The effect of on-site energies on the other hand remains, re- 
gardless of the system. We demonstrate this in part (b) of Figs. 
[TO]and[TT] which show the detrimental effect of Hi on entan- 
gled systems (types (i) and (ii) respectively). Similar to Fig. 
|6] the loss in EoF scales as an exponential in N with Gaussian 
damping in the noise parameters. Again, the influence of Hi 
has been averaged over 200 random realization using random 
numbers €i from a flat distribution, such that e; ~ eJmaxdi 
and 0<d,<l. 

Figs. [TO](b) and[TT](b) also show that for e < 0.1, EoF close 
to unity can still be achieved for all chain lengths considered. 
However, we note that for larger values of E, chains with type 
(ii) entanglement suffer significantly more entanglement loss 
than those of type (i) which are already entangled at i = 0. 
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FIG. 9: Influence of H2 (Eq.([5]l) with 7 — 0.05 on the entanglement 
of formation between the end spins of a linear 8 spin chain with input 
state JT/;™) = l/2(|0i08) + IliOg) + lOilg) + llilg)) ® IO2 • ■ • O7) 
vs. rescaled time t/ts. The first peak at i = 0.5ts = Im is 0.9996. 




12 14 



FIG. 10: For chains with type (i) entanglement (spins 1 and 2 are 
entangled at f = 0): 

(a) EoF of qubits 1 and 2 at t = f s vs. chain length A^ for three values 
of A, as labelled. Fits are according to e~^ ' " with Ao = 0.31. 

(b) Average EoF of qubits 1 and 2 a\.t = ts vs. chain length A'' for 
three values of e, as labelled. Fits are according to e ' ° with 
eo = 3.06. 



As with unentangled states, next-nearest neighbor inter- 
action is a relevant issue for entangled states. In Fig. |5] 
we showed that a 6-spin chain suffered serious fidelity loss 
(~ 20%) for relatively small A = 0.01 after about 4 periods, 
but still reached a fidelity of 0.9966 at t^i; similarly we see 
in Fig. [TT]that for the same value of A, a 6-spin chain with 
case (ii) entanglement would reach EoF of 0.99 at f = tjv/ and 
performs thus equally well. EoF of almost unity persists for 
all N considered and small A. More than 90% of the EoF is 
maintained even for A as large as 0.05 for long chains with 
type (ii) entanglement while chains with type (i) suffer signif- 
icantly more and long chains lose over 30% of their EoF for 
the same value of A. We note that for values of A > 0.05, 
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FIG. 1 1 : For chains with type (ii) entanglement (spins 1 and N be- 
come maximally entangled at t = Im): 

(a) EoF of qubits 1 and A^ at i = f a/ vs. chain length TV for 
three values of A, as labelled. Fits are according to e"'^ ' " with 
Ao = 0.63. 

(b) Average EoF of qubits I and A'^ at i = f a/ vs. chain length N 

_ jw 2 / 2 
for three values of e, as labelled. Fits are according to e ' ''o ^ith 

eo = 2.01. 



in chains with type (ii) entanglement, whilst EoF is fairly well 
maintained alt = Im, its subsequent periodicity is completely 
lost for any chain length (Fig. [12] for a demonstration of the 
same effect in an unentangled chain see Fig. [T]). 

As noted above the effect of Hj, on entangled states is dif- 
ferent from case (i) to case (ii). In Fig. [T3]we see that for 
A = 0.01 a state that is initially already entangled does not 
suffer very much and retains an EoF of over 90% at i » 4i5, 
after 4 periods, whereas a chain where the entanglement is 
created through the system dynamics suffers a loss in EoF of 
about 15% after 4 periods, at t « A.hts, as shown in Fig. [14] 

Overall, we observe that the perturbative influence of next- 
nearest neighbor interaction (Eq. (O) is the main cause for 
loss in fidelity in both unentangled states as well as states 
whose entanglement results only from their dynamics and 
which are initially unentangled, although unentangled chains 
suffer slightly more. Entangled chains with an initially en- 
tangled state on the other hand have been shown to be more 
robust against this type of defect. 

Despite these variations, our study clearly demonstrates 
that next-nearest neighbor interactions are the most damaging 
form of perturbation overall, as seen in state transfer fidelity 
or EoF. The reason for this is that the next-nearest neighbor 
interaction is the only fabrication defect or limitation in the 
set (a)-(d) that effectively opens up new channels for the sys- 
tem dynamics. The Hamiltonian H^ of Eq. (O connects chain 
sites which would otherwise be disconnected (at the same or- 
der in perturbation). It therefore facilitates a more efficient (in 
a detrimental sense) spread of excitations. The general con- 
sequence of this is more damage to transfer fidelity or EoF, 
when compared to defects (a)-(c) with the same level of noise. 
Numerical simulations (not shown) support this explanation. 




FIG. 12: Spin chains with input state \ipi„) = l/2(10iOAr) + 
lliOjv) + lOiliv) + lliliv)) ® IO2 ■ ■ -Oiv-i): 

(a) EoF of\^pi„) vs rescaled time t/ts for A = 0.05: the periodicity 
of the entanglement of formation is completely lost after the second 
peak at i = l.bts- 

(b) EoF of \tpin} vs rescaled time t/ts for A — 0.1: the periodic- 
ity of the entanglement of formation is completely lost after the first 
peak at i = 0.5ts- 



as addition of new perturbative channels for the dynamics by 
hand (as opposed to via Hs) into the full Hamiltonian can 
lead to similar results to those in, for example. Fig. [T2] It 
has been shown in Ref. ||32ll that opening new channels even 
beyond next-nearest neighbor interaction can be compensated 
for if local control within the chain is possible by adjusting the 
nearest neighbor coupling. This degree of local control may 
not always be available and/or might not be desirable, so in 
our work we consider the longer range interactions as a po- 
tential perturbation on nearest neighbor systems designed to 
produce PST. 



in. EFFECT OF NON-SYNCHRONOUS AND IMPERFECT 
INPUT OPERATIONS 
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FIG. 13: Influence of ^3 (Eq. ([6)) with A = 0.01 on the entan- 
glement of formation between the end spins of an 8 spin chain with 
input state \Tpin) = 1/72(110000000) + lOlOOOOOO)) vs. rescaled 
time t/ts- The first peak at t = 0.5is = tM is 0.9986. 
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FIG. 14: Influence of H3 (Eq. ^) with A = 0.01 on the entan- 
glement of formation between the end spins of a linear 8 spin chain 
with input state \i>,n) = l/2(|0i08) + jliOs) + lOils) + llils)) ® 
O2 • ■ • O7) vs. rescaled time t/ts- The first peak at t = 0.5is = tAf 
is 0.9976. 



tion.) The timing errors we consider are with respect to this 
reference clock. 

We will first consider the case of unentangled input states 
|0{i},n) ^nd then analyze the effects on entangled input states. 



A. Unentangled states: non-synchronous injection 



In this section we will discuss the effect of imperfect ex- 
citation injection into the device when multiple excitation 
states are considered. A possible device configuration for 
input/output of multiple excitations is sketched in Fig. [15] 
where each site in the chain is associated with a register which 
can act as input/output device. We assume that there exists a 
clock to which the machinery at both ends of the chain have 
reference. (Without such a clock even simple PST could not 
operate, as extraction has to be timed with respect to injec- 



Let us consider the injection and transfer of a state |0{i},„). 
An important question to ask is how and to what extent a 
time delay in input operations would alter the transmission fi- 
delity. This question is particularly important for multiple ex- 
citations: for a single excitation a delayed injection would not 
alter the overall state evolution, but there would be merely an 
overall time shift. For multiple excitations however, it may oc- 
cur that during the preparation of a state \(f>{iyn)^ with n > 1, 
not all input sites are accessed at exactly the same time. As 
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FIG. 15; sketch of possible multiple qubit input/output device. In 
the example each spin correspond to a specific register, while the 
active input registers are the i and j registers and the information is 
transferred to the N ~ i -{-1 and N — j -\-l via the mirroring rule. 



a consequence a system which is supposed to be prepared in 
a two-excitation state may start evolving as a one excitation 
system if the two required excitations are not injected in syn- 
chrony. 

Let us focus on the latter case. When considering such a de- 
layed input there is a finite probability that the second spin k 
we want to inject an excitation into is already occupied. The 
result of this scenario is dependent on the injection mecha- 
nism, so we will consider the two main possibilities: (I) spin 
excitation via a Rabi-flopping control pulse on spin k - appli- 
cable for example to systems in which excitations correspond 
to ground state excitons confined within a quantum dot chain 
or to flipping the spin of an electron already confined within 
the chain - and (II) injection via SWAP operation or injection 
of an additional particle in the chain. The latter may corre- 
spond e.g. to the scenario in which the state in the qubit of reg- 
ister k closest to the spin chain (see Fig. [TSl l is swapped with 
the state in the chain site k - e.g. via a train of laser pulses 
in the case of exciton qubits (see II81I ) - or to the scenario 
in which the main computation occurs via coherent electron 
transport in quantum wires (such as in ||33ll ), each connected 
to a spin chain site. 

In case (I) injecting an excitation corresponds to apply- 
ing a TT-pulse using the Hamiltonian Hj^ ~ Vl{t)\0){l\k + 
ri*(i)|l)(0|fc with Q. the Rabi frequency (we assume that 
qubits can be manipulated on an individual basis). Accord- 
ingly, at the delayed injection of the second excitation, the 
system state will evolve as 



^C{^,}.l\^{^,}.l) 



3 

with d 






,}.2\<t>{t,.M},2)+do\(f>{}.o) 



{ij,i^},2 = C{j,},i for j ^ k and do 



(9) 
_C{i,}.i- Here we 



have explicitly displayed the set of indices {i}. The last term 
in Eq. ^ corresponds to the error induced by having a non- 
zero probabiUty |c{ij.} ip of an excitation present in spin k, 
and translates, after injection, in the probability of having no- 
excitations at all present in the chain (see Fig. [T6la)). If the 
desired dynamics is such that a certain site j has unit prob- 
ability to contain an excitation at a later time t though, the 
system can be refocused by measuring site j at i: a result of 
"excitation present" would collapse the system state into two- 
excitation dynamics, which is the closer to the desired dynam- 



ics the smaller the delay; a result of "no-excitation present" 
would imply that the chain contains indeed no excitation at 
all. We underline that the latter result could be used as a pro- 
tocol to reinitialize the chain itself. 

In case (II) the scenario is very different. Should the first 
excitation already have a non-zero probability of occupying 
site fc, the second excitation will have a related probability of 
remaining in the register. The register would then become en- 
tangled with the spin chain. The injection process can in fact 
be described as follows (where we assume that the presence 
of the excitation in site k is the only cause of failed injection) 

2^C{i^}^l|0{i^}4) (g) \l)register ^ 
j 

C{*fc},ll<^{jfc},l) ® ^register, (10) 



{i,,*.},2=qi,}.iforJ7^A:. 



with d 

The simplest way to destroy this unwanted entanglement is 
to measure, after injection, whether the second excitation is 
still present in the register or wire. By this measure, we re- 
move the entanglement between register and device, but we 
also get to know exactly what state the device itself is in: if 
the measure outcome is "no excitation in the register" the in- 
jection has been successful and the spin chain now follows 
a two excitation evolution which is the closer to the desired 
dynamics the smaller the injection delay has been. This is de- 
scribed in Fig. [T6lb). If the outcome of the measurement is 
"excitation in the register" the chain will continue to evolve as 
a one excitation system. However in the latter case we have 
also re-set the chain ready for trying again the injection of the 
second excitation at the earliest convenient time. 

In Fig. [TtJ we give an overview of various possible delay 
scenarios, assuming input by SWAP operation with error cor- 
rection by measurement of the injecting register or wire per- 
formed immediately after the attempted injection of the sec- 
ond excitation. In panel (a), the delay of the second excitation 
is equal to 0.15^^/. Even with the error correction measure- 
ment we see an impact on the system as the delay between the 
two injections puts a cap on the fidelity of the revivals of the 
desired two-excitation input state. The larger the delay (up to 
the mirror time tj\/), the more serious the impact. Frame (b) 
demonstrates that even if timely simultaneous injection were 
not possible, the ability to inject the second excitation an in- 
teger number of periods ts after the first one would allow for 
perfect revivals. In contrast, the third frame (c) shows that a 
second injection exactly Im away from the ideal case would 
result in the total decay of the desired state. However, we have 
to keep in mind that at t == tjv/ where \(t>{i}.i) h^s the lowest 
fidelity, the corresponding twin state will be achieving perfect 
fidelity and that therefore a new state will emerge, we illus- 
trate this in Fig. [18] at a time equal to an odd multiple of 
tjv/, the initial input state 1 100000) has zero fidelity, while its 
twin state jOOOOOl) has unit fidelity. Injection of the second 
excitation into the second site (which would have resulted in 
the ideal state 1 110000) if there had been no delay) results in 
lOOOOOl) -^ lOlOOOl). This state then continues to evolve 
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FIG. 16: Fidelity with input of the second excitation in a 6-spin chain 
with O.Um delay vs. rescaled time t/ts'- 

(a) shows injection by Rabi-flopping, where the resulting error re- 
mains in the system as the zero vector 1 000000). At t — ts the 
state of the first spin is measured and an excitation is found. This 
allows for refocusing the system in the two excitation subspace only. 
The peak of 1 110000) at multiples of ts plus delay is 0.7807. If at 
t = ts the second spin is measured instead, the peak fidelity be- 
comes 0.7203. 

(b) shows injection by SWAP operation where the error remains in 
the system in the one excitation subspace. At i = is the state of the 
register is measured and no excitation is found. This projects the sys- 
tem in the two excitation subspace only and disentangles the chain 
from the register. The peak of 1 110000) at multiples of ts plus delay 
is 0.9870. 
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FIG. 17: Effects of delayed input of a second excitation on a 6-spin 
chain with SWAP operation given by fidelity vs. rescaled time t/ts- 

(a) input delayed by 0.15iA/. 

(b) input delayed by an integer number of ts (here: Its). 

(c) input delayed by an odd multiple of f Af (here: Um)- 

For (a) the maximum recumng fidelity of |000011) is 0.9313, at odd 
integer multiples of tm plus the delay. For (b), due to the complete 
period delay, jOOOOll) emerges with unit fidelity, whereas for (c) the 
occurrence is negligible. 



with its twin state 1 100010), both of them alternately reaching 
perfect fidehty. The additional features on the plots of vectors 
with two excitations, as well as the narrowing of the main peak 
are due to excitation of the particular superpositions of energy 
eigenstates in the two-excitation subspace that correspond to 
our initial state. This shows how we can use delayed input to 
transfer states which are different from the "twin" sites of the 
input sites while still assuring PST. This may be useful if e.g. 
not all sites are accessible for input operations. 
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B. Unentangled states: Imperfect input operations 



Imperfect injection can occur even without delays between 
different excitations' injections due to other device imperfec- 
tions. A typical example could be the possibility that in the 
SWAP injection scenario one of the two excitations is par- 
tially reflected into the wire. This may occur in a device in 
which the spin chain is formed by gate-defined quantum dots 
where electrons are injected via (computational) wires, with 
the reflection being caused by imperfect lowering of a poten- 
tial barrier between the wire and the selected spin chain site 



FIG. 18: Fidelity after second input to a 6-spin chain is delayed by 
an odd multiple of ti\/ (here: Sti/) vs. rescaled time t/ts- The peak 
of lOlOOOl) at t = 2.5fs reaches unity. 



k- In this case the injection would be described by 

I0{},O/ ® |l/registerfc ® \^)registerj ~^ 



d{ij,ik},2\'t'{ij,ik}a) ® \^)registerk ® \Q)registerj 



C{^^},M{i,}, 



ii; 



register]^ 



®|0) 



register j i 



(11) 



where the first spin has been perfectly injected at site j, while 
there is a reflection probability |c{i j ip associated with spin 
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FIG. 19: Maximum value of the first entanglement of formation peak 
at tM vs. input delay for an 8-spin chain with initial input state 

\i>,n) = l/2(|0i08) + IliOg) + lOilg) + llilg)) ® |02 • • • O7) for 
both Rabi flopping and SWAP operation type of injections. 



FIG. 20: Maximum value of the first entanglement of formation peak 
vs. input delay for a 6-spin chain with initial input state {il^in) ~ 
l/^/2(| 100000) + lOlOOOO)) for both SWAP operation and Rabi 
flopping type injections. 



k. In this scenario measuring the absence (presence) of excita- 
tion in the register after injection would ensure that the chain 
is undergoing exactly the desired dynamics (or that the chain 
is ready for re-injecting the second excitation at the closest 
suitable time). The above scenario can be straightforwardly 
extended to the case in which a finite probability of reflection 
is associated with both injection sites. 

While in this paper we are only discussing delays at the in- 
put stage, similar problems may arise of course at read-out, 
should the extraction of a state covering multiple sites not be 
as timely as we would hope it to be. It is also worth noting 
that any peak following delayed input is shifted forward in 
time linearly with increasing delay. 



When we consider two adjacent injection sites, there is in- 
stead a clear discrepancy between the decay in entanglement 
of formation in a chain with input via Rabi flopping or SWAP 
operation (Fig|20]i. Again, we assume that there is an instan- 
taneous SWAP correction protocol applied, while refocussing 
after the Rabi-flopping is not possible. As in the case of non- 
synchronous input in unentangled states, a delay in a SWAP 
type injection perturbs the system far less than using Rabi 
flopping, and shows virtually no loss in the amount of entan- 
glement in the system, even for delays as large as 0.05ij\/. If 
instead Rabi flopping is used, the decay in EoF is similar to 
that of Fig. [19] 



C. Entangled states 

When considering imperfect or non-synchronous input and 
entanglement transfer, we see a similar effect to that on un- 
entangled states. Again, as the overall desired system state 
is more complicated, we monitor the evolution of the amount 
of entanglement in the relevant sub-system, instead of indi- 
vidual vectors. In Fig. [19] we show the evolution of the first 
EoF peak after the delayed input has taken place, for case (ii) 
where entanglement is created from an initial product state. 
In this setting it is not possible to implement the refocusing 
protocol for the Rabi flopping, but we are assuming that the 
confirmed injection protocol for the SWAP operation is imple- 
mented instantaneously, just after the injection attempt. Due 
to the fact that the two injection sites in the state described 
in this figure are at opposite chain ends, there is virtually no 
difference between injection by Rabi flopping and injection 
by SWAP operation. As long as the second injection is only 
delayed by a small fraction of tnj, the first excitation will not 
have spread out far enough yet to affect the second injection 
site. We see that even for delays as large as O.Um, the system 
retains over 90% of its possible entanglement. 



IV. CONCLUSIONS 

In this paper we have considered a variety of physically rel- 
evant perturbative factors in spin chains. We have investigated 
in detail their effects on information transfer and entangle- 
ment generation and transfer For different forms of pertur- 
bation, the results of our extensive numerical studies on the 
quality of transfer can be captured with a straightforward an- 
alytic expression that demonstrates exponential damping with 
the number of spins on the chain and Gaussian dependence on 
the relevant perturbation parameter. This expression provides 
a simple tool for estimating the efficiency of a chain under 
the action of perturbations. We have also related our dimen- 
sionless perturbation parameter scales to actual parameters for 
candidate spin chain realizations such as quantum dots with 
exciton qubits and graphene dots with spin qubits, providing 
a calibration of our estimator tool for these experimental sys- 
tems. 

We have considered the transport of both unentangled states 
as well as entangled states via spin chains, subject to various 
forms of perturbation, and a general conclusion is that next- 
nearest neighbor interactions are the most damaging to trans- 
fer fidelity or EoF. The reason for this can be traced back to 
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the fact that the next-nearest neighbor interaction is the only 
perturbation that effectively opens up new 'channels' for the 
system dynamics. By connecting sites which would other- 
wise be disconnected (at the same order in perturbation), it 
allows a more efficient (in a detrimental sense) 'spread' of ex- 
citations, and consequently drastically diminishes the occur- 
rence of quantum coherent effects such as revivals and PST. In 
particular, we have seen that the introduction of next-nearest 
neighbor coupling may lead to non-negligible quality loss af- 
ter only a few periods. Systems with simpler input states 
are generally slightly less affected than those of a slightly 
more complex nature but as different perturbation factors af- 
fect different excitation subspaces, there is no clear advantage 
of one particular type of state in terms of robustness. How- 
ever, for all fabrication defects considered we have found that 
the transport during the first period remains of high quality 
for perturbation amplitudes of the order of few percent, while 
the periodicity itself may be destroyed for next-nearest neigh- 
bor couplings above about 5%. This is of particular concern 
for schemes based on dipole-dipole interactions, where, for 
roughly equidistant sites, the next-nearest neighbor interac- 
tions are of the order of 10%. Nevertheless, due to the Gaus- 
sian dependence of transport efficiency on perturbation ampli- 
tudes, provided that the perturbations are kept below these few 
percent thresholds, spin chains are demonstrated to be very 
good candidates for the implementation of solid state quantum 
information processing devices, robust against various forms 
of perturbation. 

Additionally, we have seen that the effect of imperfect in- 
put or injection operations leads to a permanent (but constant) 



loss in information transport quality and entanglement gener- 
ation. For input delayed by a time approaching the mirroring 
time tm the transport of the intended state is replaced by a 
new input state, which is then subject to the system dynam- 
ics as per usual. Furthermore, there is a fundamental differ- 
ence between the two input methods we have considered. If 
Rabi flopping is used and the injection of the second excita- 
tion is delayed, the error induced by the imperfect input re- 
mains in the system, except if a site is expected to be in state 
|1) at a known time, when measurement of this site can re- 
cover dynamics close to the ideal case. In the case of input 
by SWAP operation, the error can be dramatically reduced by 
subsequent measuring of the environment. When considering 
SWAP type injection, unentangled states and type (i) entan- 
glement states are remarkably robust against non-synchronous 
injection, with hardly any loss in fidelity even for large delays. 
Type (ii) dynamically-generated entangled states are instead 
more affected due to their more complex set-up. Even so the 
loss in EoF remains less than 10% for delay values up to 10% 
of the mirroring time tnj- 

Our studies demonstrate quantitatively the criteria that need 
to be met, in terms of perturbation scales and injection errors, 
for imperfect spin chains to work efficiently as quantum infor- 
mation transfer and entanglement transfer/generation devices. 
With modest errors at or below the few percent level, spin 
chains prove to be good and robust devices, which is very en- 
couraging for future experiments on these systems. 
RR was supported by EPSRC-GB and HP IDA acknowledges 
partial support by HP. IDA and RR acknowledge the kind hos- 
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